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1 Introduction 



Flavor changing neutral currents (FCNC) and charge parity (CP) violat- 
ing phenomena are some of the best probes for physics beyond the Stan- 
dard Model (SM). All existing measurements so far are consistent with the 
SM picture involving the Cabibbo-Kobayash-Maskawa (CKM) matrix as 
the only source of flavor violation. In the SM, FCNC are absent at tree 
level, appear at one loop, but they are effectively suppressed by the GIM 
mechanism and small CKM angles. In supersymmetry, there is no similar 
mechanism to suppress the loop contributions to either flavor or CP vio- 
lating phenomena. Experimental studies of flavor physics, especially in B 
decays, appear essential for the understanding of the mechanism for super- 
symmetry breaking. With increased statistical power of experiments at B 
factories, rare B decays will be measured very precisely. 

In this paper we investigate the relevance of new physics in the semilep- 
tonic inclusive decay B — > X s l + l~ in a fully left-right supersymmetric 
model. Investigation of the process b — > 57 in this model has shown dis- 
tinctive signs from the minimal supersymmetric standard model (MSSM) 
scenario ||T|. An analysis of B — > X s l + l~ would provide some complemen- 
tary information. The semileptonic decay B — > X s l + l~ is a benchmark 
of charmless 6-decays with strange particles in the final state. The pro- 
cess is experimentally clean, but the expected SM branching ratio in the 
— 1CT 5 region makes it not easily detectable at B-factories. There- 
fore it provides for excellent opportunities to test physics beyond the SM. 
It also offers more detailed information about the flavor structure of the 
model, it provides a good test of the structure of the Zbs vertex, making 
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it particularly well suited to distinguish the left-right symmetry over the 
MSSM. 

Experimentally, BELLE has recently announced the first evidence for 
the exclusive process B — > K*l + l~ . BABAR and BELLE have upper 
bounds for B — ► (K,K*) + (e + e~ , fi + fi~) which are very close to the SM 
estimates. Experimentally, the bounds on the branching ratios are [fj 

BR(B -> X s /i+/0 < 19.1 x 1(T 6 @ 90% C.L., 
BR(B -> X s e + e-) < 10.1 x 10~ 6 @ 90% C.L., 

typically a factor of 3 away from the SM estimates, where the next-to-next- 
to leading logarithmic (NNLO) calculations to 0(1/ oc s ) || have appeared 
recently 

BR(B -> X s n + iT) < (4.15 ± 0.70) x 10" 6 , 
BR(B -> X s e + e") < (6.89 ± 1.01) x 10~ 6 . 

Here we restrict ourselves to the analysis of inclusive processes only: in 
the case of exclusive decay rates, hadronic matrix element uncertainties 
obscure model predictions. 

Semileptonic charmless B decays have been studied previously by many 
authors in the framework of supersymmetric models with a universal soft 
supersymmetry breaking terms [|J. Recently an analysis of SUSY models 
with non-universal soft breaking terms at the grand unification scale has 
appeared in Ref. 0. Although attempts have been made to reconcile 
b — ► S7 with right-handed 6-quark decays 0, a complete analysis of B — > 
XJ + l~ for a fully left-right supersymmetric model is still lacking. 

The Left-Right Supersymmetric (LRSUSY) models [0, ||, based on the 
symmetry group SU{2) L x SU{2) R x U(1)b-l, incorporate the advantages 



of supersymmetry with a natural framework for allowing neutrino masses 
through the seesaw mechanism |J. Various other scenarios incorporate 
some forms of the left-right symmetry within supersymmetry. LRSUSY 
models have the attractive feature that they can be embedded in a super- 
symmetric grand unified theory such as SO(10) [p70| , while not bound by 
lepton quark unification. They would also appear in model building realistic 
brane worlds from Type I strings. This involves the left-right supersymme- 
try, with supersymmetry broken either at the string scale Msusy ~ 10 10 ~ 12 
GeV, or at Msusy ~ 1 TeV, the difference having implications for the gauge 
unification [fllll . 

In this paper we study all contributions of the LRSUSY model to the 
branching ratio and the asymmetry of B — > X s l + l~ at one-loop level. 
The decay can be mediated by the left-handed and right-handed W and Z 
bosons, and by charged Higgs bosons as in the nonsupersymmetric case, but 
also by charginos, neutralinos and gluinos. The structure of the LRSUSY 
model provides a significant contributions to the decay B — > X s l + l~ from 
the right-handed squarks and an enlarged gaugino-Higgsino sector with 
right-handed couplings, which is not as constrained as the right-handed 
gauge sector in the LRSUSY model. We anticipate that these could con- 
tribute a large enhancement of the decay rate and would constrain some of 
the parameters of the model. 

The paper is organized as follows. We describe the structure of the 
model in Sec. II, with particular emphasis on the gaugino-Higgsino and 
squark structure. In Sec. Ill, we give the supersymmetric contributions in 
the LRSUSY model to the decay B — > X s l + l~ . We confront the calcula- 
tion with experimental results in Sec. IV, where we present the numerical 
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analysis to constrain the parameters of the model for two scenarios: one 
with the CKM flavor mixing only, the other including supersymmetric soft 
breaking flavor violation terms. We reach our conclusions in Sec. V. 

2 The Model 

The LRSUSY electroweak symmetry group, SU(2) L x SU(2) R x U{\)b-l, 
has matter doublets for both left- and right-handed fermions and their cor- 
responding left- and right-handed scalar partners (sleptons and squarks) || . 
In the gauge sector, corresponding to SU(2)l and SU(2) R , there are triplet 
gauge bosons (W + 1 W~ ,W°) L , (W + , W~, W°) R , respectively, and a sin- 
glet gauge boson V corresponding to U(1)b-l, together with their su- 
perpartners. The Higgs sector of this model consists of two Higgs bi- 
doublets, |, 0) and |, 0), which are required to give masses to 

the up and down quarks. The spontaneous symmetry breaking of the group 
SU{2) R x U(1)b-l to the hypercharge symmetry group U(l)y is accom- 
plished by giving vacuum expectation values to a pair of Higgs triplet fields 
A^(l,0,2) and A R (0, 1,2), which transform as the adjoint representation 
of SU (2) R . The choice of two triplets (versus four doublets) is preferred be- 
cause with this choice a large Majorana mass can be generated (through the 
see-saw mechanism) for the right-handed neutrino and a small one for the 
left-handed neutrino ||. In addition to the triplets iR , the model must 
contain two additional triplets, 5^(1,0,-2) and <5r(0, 1,— 2), with quan- 
tum number B — L = —2, to insure cancellation of the anomalies which 
would otherwise occur in the fermionic sector. The superpotential for the 
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LRSUSY model is 

W LRSUSY = \ifQ T r 2 ^ t r 2 Q c + h?L T r 2 ^ l T 2 L c + i{h LR L T r 2 A L L 
+ h LR L cT r 2 A R L c ) + M LR [Tr(A L 5 L + A R 5 R )} 
+ Hi j Tr(T 2 <S>jT 2 $ j ) + W NR , 

where W^ R denotes (possible) non-renormalizable terms arising from higher 



scale physics or Planck scale effects [12]. The presence of these terms insures 



that, when the SUSY breaking scale is above Mw R , the ground state is Im- 
parity conserving ||13|| . 

The neutral Higgs fields acquire non-zero vacuum expectation values 
(VEV's) through spontaneous symmetry breaking 

< a> « - ( vi o ) • - <•>■* - ( r <> ) ■ 

(&)u,d cause the mixing of Wl and W R bosons with CP-violating phase u, 
which is set to zero in the analysis. The non-zero Higgs VEV's break both 
parity and SU(2) R . In the first stage of breaking, the right-handed gauge 
bosons, W R and Z R acquire masses proportional to v R and become much 
heavier than the SM (left-handed) gauge bosons Wl and Zl, which pick 
up masses proportional to n u and Kd at the second stage of breaking. 

In the supersymmetric sector of the model there are six singly-charged 
charginos, corresponding to X L , \ R , (j) u , 4>d, Af, and A^. The model also 



has eleven neutralinos, corresponding to Xz, \z', Ay, 4>1i, 0°2> 4>dii 0°2> 
A° L , A° R , S L , and 5%. Although A^ is not necessary for symmetry break- 
ing [TJ]] , and is introduced only for preserving the left-right symmetry, both 



A L (A L ) and its right-handed counterparts A^ (A^ ) play very impor- 
tant roles in lepton phenomenology of the LRSUSY model. The doubly 
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charged Higgs and Higgsinos do not affect quark phenomenology, but the 
neutral and singly charged components do, through mixings in the chargino 
and neutralino mass matrices. We include only the A# contribution in the 
numerical analysis. 

The supersymmetric sources of flavor violation in the LRSUSY model 
come from either the Yukawa potential or the trilinear scalar couplings. 
The interactions of fermions with scalar (Higgs) fields have the following 
form 

C Y = KQl^uQr + h d Q L <£> d Q R + h u L L <Z> u L R + h e L L <Z> d L R + H.c, 
C M = th LR (L T L C- 1 T 2 A L L L + L T R C- 1 r 2 A R L R )+ H.c., 

where h u , h d , h u and h e are the Yukawa couplings for the up and down 
quarks and neutrino and electron, respectively, and \il R is the coupling for 
the triplet Higgs bosons. The left-right symmetry requires all h-matrices 
to be Hermitean in the generation space and \il R matrix to be symmetric. 
In the universal case, there is no intergenerational mixings for squarks and 
the only source of flavor mixing comes from the CKM matrix. We will 
analyze this case first. Next we will look at the case in which intergener- 
ational mixings in the squark sector are permitted and consider the effect 
of intergenerational mixings on the rate of the process B — ► X s l + l~ . 

3 The analytic formulas 

The effective Hamiltonian for the decay B — > X s l + l~ at the scale fx in the 
LRSUSY model is given by 

AC 



The operators relevant to the process b — > s/ + Z in the LRSUSY model are 

Q 7 = m b (q u /q 2 )(sia^P R b)(h^), 

Q 7 = m b (q u /q 2 )(sia^P L b)(l'y fl l), 

Q9 = {s^P L b){l^l), 

Q9 = (57^6) ^Z), 

Q10 = (s7 At P L 6)(l7 M 7 5 /), 

Q10 = (s^nPR^ih^l)- 

The Wilson coefficients Cj and Ci are initially evaluated at the electroweak 
or soft supersymmetry breaking scale, then evolved down to the scale \x. In 
the SM and constrained SUSY models, Qi contributions are generally sup- 
pressed by 0{m s /m^) compared with the contributions from Qi. However 
this is not the case in generic SUSY models such as non-universal models. 
In Ref. H the operator Q 7 was included in the analysis. Due to the left- 
right symmetry, we are motivated to consider all contributions from both 
chirality operators. 

The decay b — > sl + l~ can be mediated by either the photon or the 
Zl, Zr bosons, or it can proceed through the box diagrams. As in the 
MSSM, the Zl boson contributions dominate where there is explicit SU {2)l 
symmetry breaking, and the Zr boson contributions are important where 
there is explicit SU(2)r symmetry breaking, i.e., both cases in which left 
and right squarks occur in the same loop diagram. In these cases, the Z 
diagrams are enhanced by m?/M| with respect to the photon graphs. This 
could be an order of magnitude for the regular Zl boson, but only of order 
1 for the Z R boson. We describe these contributions in detail below. 
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Figure 1: Penguin diagrams which induce the decay b — > s/ + Z in the LRSUSY model. 
The outgoing photon and Z boson lines can be attached in all possible ways. 

3.1 The photon monopole and penguin graphs 

The penguin monopole and dipole graphs are presented in Fig. [I]. We first 
give the contributions for the constrained model, i.e., assuming the CKM 
matrix is the only source of flavor violation. The left-handed monopole 
contributions are given by 
a-^oi 1 
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2 M^ L 



K* ts K tb {xtw[QuF 7 (x t w) + Fg(xtw)] 



2 Mx^ 

3 x t w - 1 
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LL 



Kl s K tb x t H cot 2 P[Q u F 5 (x t H) ~ F 6 (x tH )}, 



A 



LL 



OL-^a 1 

5 6 

Z ^t( G ul - h ur)( g ul - Hi 



j=l k=l m u k 



r*jks\ 

l UR ) 



x[Q u F 6 (x~ 7Uk ) - F 5 (x 



A L 



LL 



LL 



A x° 



1/1 flj 

3=1 k=l d k 



The monopole contributions with Wl replaced by Wr, as in the left-right 

M 2 

symmetric model (LRM), are suppressed by M ^ L and thus negligible. The 



(7) 
(8) 

(9) 
(10) 



-2a s aC{R)Q d £ — r^r^^-jj, 
fc=i m d k 

-awaQ^Y: ^-(v^Gjg z - ^(^Cg* - ^)^(^o 4 ). (11) 
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only RR monopole contributions then come from the supersymmetric sector 

a rr = -%«EEi(GS-C)(G^ s -^ s ) 

3=1 k=l m u k 

x[^(^ 7 , fc )-F 5 (^-, fc )], (12) 

Af R = -2a a aC(R)Q d j2^T% R r* R F 9 ( X9 ri, (13) 

k=i m d k 

A R o R = -a w aQ d £ £ ^(^^ - ^(v^GgS - flSS)F 6 (^). (14) 

j=l k=l m d k 

The dipole LR and RL contributions can be obtained by multiplying the 
contributions to the decay b — > 57 by the factor V^ra 

^sm = °T ^ K* ts K tb 3x tw [Q u Fi(x tw ) + F 2 (x w )], (15) 
2 M w L 

A L H R - = aa Q W * K* ts K tb x tH {cot 2 f3[Q u Fi{x tH ) + F 2 {x tH )] 
1 M w L 

+ [Q u F 3 (x tH ) + F 4 (x tH )]}, (16) 

6 1 m- 
A L g R = -2aa s Q d C(R) ]T — r {T'S L T^ s L F 2 (x §(Ik ) - -^T>g R r*£ s L F 4 (x- gd ~ k )}, (17) 

k=i m d k mb 

j=l k=l m u k 3 J 

m x7 



+^r( G UR - Hh kb L )(GfL - H*J k R s )[F,{x^ k ) + Q u F 4 (x- T , k )]}, (18) 



3=1 k=l m d k 



+^(x/2G& - HjgJiVxTjSl - H*^)F 4 (x~o dk )}, (19) 



and 



6 1 777 ~ 

Af = -2aa s Q d C(R)^^{r% R ri k R F 2 (x~ d ^ (20) 

k=i m d k nib 

Af = "««^ E E - - ^ S )[fl(^:,J + Q U F 2 (X.,J] 



3=1 k=l m u k 



m x 



^( G ul ~ Hb kb R )(Gt R s ~ H*jtm(x~-J + Q„F 4 (x X7fifc )]}, (21) 
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j=i k =i m d k 

where vertex mixing matrices G, H, G and if are defined in the Appendix. 
The convention x a b = m-a/ml * s used. C(R) = 4/3 is the quadratic Casimir 
operator of the fundamental representation of SU(3)c- 

3.2 The Zl, Zr penguin graphs 

The process b — > sl + l~ is also induced by the effective couplings of the Zl 
and Z R . The diagrams are analogous to the photon graphs where the pho- 
ton is replaced by the Z propagator. The amplitudes for the Zl mediated 
graphs are 

a 2 1 

A Z s \l = cos2 dw KK tb x tw F g (x tw ), 

a 2 1 

A h L - = ITTTi ^z- K* ts K tb cot 2 (3 XtwXtHiF^Xtn) + F A (x tH )], 

4 IVl Zl COS u\v 

m Zl cos v w h k=1 



m=l 



6 5 



aZ l CHyyCH 1 \ -< \ ^ /^yjkb Tjjkb\ / ^y*ihs u*ihs\ 

K m=l 

xUaU^-Goix^x^V^]}, 
z iw ZL cos t/ w hjk=li:j=1 

K m=l 

x(N t4 N* 4 - N l5 N* 5 ) - Go(x^ k ,x^ k )(N* 4 N j4 - N* 5 N j5 )]} , 
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where we have included in the expressions both cases in which one vertex is 
gaugino and the other Higgsino, and the case in which we have two gaugino 
vertices. Similarly we obtain, for the Zr mediated graphs 



W 

? / 1 DR L DR 

m=l 

6 5 



A z * = -a w a s C(R) J £ T dr T dr E V^T^G^ 

M'Z„ COS our , , - , , 



S> X d h g)' 



A Z R ay/® COS 19\y x - frtjkb TTjkb\ (r i*ihs zj*ihs\ 

A X~ ~ 9~ A,f2 pnQ 2/J 2s UR~ ^ULJV^UR ~^UL ) 

2 M Zr cos t>w htk=liJ=1 



3 



I m=l 

xU^U; 2 -G (x x--JV^V j2 \ 
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aZ r .... ayj/Q cos 26^ x - ST (^/ori kb W kb \(^f 7 >r* ihs n* ihs \ 
2 M z cos «w fcfc=li =1 
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As in the MSSM [|J, when the Zl,r bosons are exchanged, gluino- and 
neutralino-induced contributions to the total amplitudes are suppressed by 
0(q 2 /M Zl ) with respect to the photon penguin contributions. In addi- 
tion, the Zfj contribution is suppressed with respect to the left-handed one 
by <D{M%JMI R ). 

3.3 The box diagrams 

The box graphs are presented in Fig. |^. The explicit contributions are 
given by 

A a SM = -^j^K* s K tb [G(x m ,0)-G(0,0)}, 

n 2 6 5 1 
aLD _ u ff \ ~» \ ^ x t/~tjkb Tjjkb\f ^*iks Tj*iks\ 

A X~ ~ ~T~ 2^ Zs —Q-V^UL- ^URtV^UL - n UR ) 

* h,k=l i,j=l " l x~ 
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Figure 2: Box diagrams contributing to the decay b — > s/ + Z in the LRSUSY model. 
Clashing arrows on the fermion lines indicate a Majorana mass insertion. 
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i*ihl rijhl 



*iks \ 
ODRJ 



^OER^OER 



XG '( X d k xV X hxV X X°X°^ ~ 2G 0ER G 0ER X /W^( x d kx 0, X l h xV X m) 



(34) 



The right handed supersymmetric contribution is 
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A x° 
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E E 



w 



4 ," 1 mt- 

h,k=li,3=l x 
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UR n ULJ\ u UR n VL J 



*iks\ 



xG NR G NR G ( X u k x7 ,X i>hX7 ,X x7x7^ 
2 9 6 | 



^OER^OER 



yE E ^(V2G^ R -Ht D b L )(V2G^-H^l) 
xG "( x 4x°> x kx - x x°x°) ~ 2G lf R G i h ER^x^F\ x d k xO, x x°x?) 



2 9 6 



1 



(35) 



(36) 



, E E 7^(^ G 0DR- H 0Dl)(^ G *0DR- H t 
Z i,j=l k,h=i rn x° 



*iks \ 
ODLj 



^OEL^OEL 



xG"(xj^o, a^o, ^°x°) - 2G* OEL G J OEL Jx^.F'(x^o, x !h ^,x^o) 



• (37) 



All the relevant vertex and loop functions are listed in the Appendix. 



13 



3.4 Branching ratios and asymmetries 

Putting all the above contributions together, we write the total amplitude 
at the M w scale as 

A tot (M w ) = C 7 (M W )Q 7 + C 7 (M W )Q 7 + C 9 (M W )Q 9 + C 9 (M W )Q 9 
+ C W (M W )Q W + C 10 (M W )Q 1Q , 

where 

A tot (M w ) = C 7 (M W )Q 7 + C 7 (M W )Q 7 + C 9 (M W )Q 9 + C 9 (M W )Q 9 
+ C 10 (M W )Q 10 + C W (M W )Q 10 , 

where 

C 7 (M W ) = Al R , 
C 7 (M W ) = Al L , 

C 9 {M W ) = A\ L + (~ + sin 2 w )A z - + ^(A a SM + A^ + A^ + A^'), 
C 9 (M W ) = Al R + l -{l + s^e w )A z - + l -{A^ + Af + Af'), 

c w (m w ) = V--i(^ M + 4 n + 4o D -4o n '), 

C W {M W ) = sm 2 e w A z « + ±(A^ + A??-A^'), 
with 

A" 1 A LL 4- 4 LL -I- A LL -A- A LL 4- A L £ 

4 7 A RR I /t -Ri? I /t -Ri? 

71 Ri? — A g "T X° ' 

A 1 A LR 4- 4 Li? 4- 4 Li? 4- 4 LjR 4- A L ? 

LR — A SM + A H~ ^ A g + ^x" + X° ' 

47 4.RL I ARL , 4.RL 

^RL y ^X" ^ ^X ' 

^ L = A§ l m + A%l+A? l +A£l+A$;, 

A Z R = A Z R+A Zn +A p_ 
14 



As the experimental results on semileptonic decays B — > X s l + l~ fit the SM 
well, new physics effects can be parameterized by Ri and Ri defined at the 
electroweak scale as 

ri r~<SM ri 
_ - ~ _ Uj 

^ ~~ qSm ' ^ ~~ qsm ■ y ^) 

Note that there are no contributions to Ci in the SM. The non-resonant 
branching ratios can be expressed in terms of the parameterization as || 

BR(B -> X s e + e") = 7.29 x l(r 6 [l + 0.35^7 + 0.179i? 9 + 0.714i? 10 

+0.0947(^ + R 2 7 ) + 0.045(^ + R 2 9 ) + 0.357(i^ + R 2 1Q ) 
-0.0313(R 7 R 9 + R 7 R 9 )}, (53) 

BR(B -> X s /i + /i") = 4.89 x 1(T 6 [1 + 0.0982^7 + 0.264^9 + 1.07i?io 

+0.0491(i^ + R 2 ) + 0.0671(^ + R 2 9 ) + 0.535(i^ + R 2 W ) 
-0.0467(^9 + #7^9)] • (54) 

If Ri and Ri are set to zero, the SM values for the semileptonic decays are 
recovered in these formulas. Resonant contributions were studied in Ref. 
|I~5| and these can be avoided by excluding some special areas from the 



phase integration regions in the dilepton invariant mass. 

We also consider the lepton pair energy asymmetry in the decay B — > 
X s l + l~ defined as 

N(E l ->E l+ )-N(E l+ >E l -) 
~ N(E t - >E l+ ) + N(E l+ >E l -Y 

where, for instance, N(Ei- > Ei+) is the number of the lepton pairs where 
the negative charged lepton is more energetic than the positive charged 
lepton in the B rest frame. The energy asymmetry is equivalent to the 
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ordinary forward-backward asymmetry. In a configuration where l + is scat- 
tered in the forward direction, kinematically, in the dilepton center-of-mass 
frame, it is implied that E t + < E t - in the B rest frame. With the above 
parameterization, the energy asymmetry is found to be 
0.48 x 1(T 6 

A i + i~ = ——— [1 - 0.625^7 + 0.884^9 + 0.911^0 

^BRltf — > A s l+1 ) 

-0.625(R 7 R W + R 7 R 10 ) + 0.884(i? 9 i? 10 + RgRia) 
-0.00882(^ + R 2 10 )\, 

where R B r(B -> X s l+r) = b§SSx^FI^ • 

4 Numerical Analysis 

We are interested in analyzing the case in which the supersymmetric part- 
ners have masses around the weak scale, so we will assume relatively light 
superpartner masses. We diagonalize the neutralino, chargino, scalar quark 
and lepton mass matrices numerically and require in all calculations that 
the masses of gluinos, charginos, neutralinos, squarks and sleptons be above 
their experimental bounds. There are some extra constraints in the non- 
super symmetric sector of the theory, requiring the FCNC Higgs boson $^ 



to be heavy, but no such constraints exist in the Higgsino sector [16|. We 



constrain the lightest Higgs boson mass to be 115 GeV [17]. 

As the first step, we assume the only source of flavor violation to come 
from the CKM matrix. This scenario is related to the minimal flavor viola- 
tion case in supergravity. This restricted possibility of flavor violation will 
set important constraints on the parameter space of the LRSUSY model. 

We then allow, in the second stage of our investigation, for new sources 
of flavor violation, coming from the soft breaking terms. In the MSSM, this 
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scenario is known as the unconstrained MSSM. We restrict all allowable LL, 
LR, RL and RR flavor mixings, assuming them to be dominated by mixings 
between the second and third squark family. 

We now proceed to discuss both these scenarios in turn. 

4.1 The constrained LRSUSY model 

By the constrained LRSUSY model, we mean the scenario in which the only 
source of flavor violation comes from the quark sector, through the CKM 
matrix, which we assume to be the same for both the left- and right-handed 
sectors (manifest left-right symmetry), as explained below. 

Before any meaningful numerical results be obtained, explicit values for 
the parameters in the model must be specified. There are many parameters 
in the model, such that it is hard, if not impossible, to get an illustrative 
presentation of calculation results. If the LRSUSY model is embedded in a 
super symmetric grand unification theory such as SO (10), there exist some 
relationships among the parameters at the unification scale Mqut- We can 
generally choose specific values for parameters at the mass scale \i = Mqut, 
then use renormalization group equations to run them down to the low 
energy scale which is relevant to phenomenology. But, for maintaining 
both simplicity and generality, we can present an analysis in which the 
LRSUSY model is not embedded into another group. Then we can choose 
all parameters as independently free variables, with the numerical results 
confronting experimental bounds directly. 

To make the results tractable, we assume all trilinear scalar couplings 
in the soft supersymmetry breaking Lagrangian as Aij = A5ij and ^ = 
/iSij, We also set a common mass parameter for all the squarks M 0UL = 
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Figure 3: BK(B — > X s l + l~) normalized to the corresponding SM values as a function of 
tan/5, obtained when m- g = 300 GeV, \i = 100 GeV, A=50 GeV, M L = M R = 500 GeV, 
moq — 300 GeV and m i = 100 GeV. The experimental constraints are also shown. 

M 0UR = M 0DL = M 0DR = m 0g . We take K% KM = K§ KM . This choice is 
conservative, and much larger values of mixing matrix elements are allowed 
in scenarios that attempt to explain the decay properties of the b quark 
as being saturated by the right-handed b ||. Our choice does not favor 
one handedness over the other, and has the added advantage that no new 
mixing angles are introduced in the quark matrices. 

We investigate first the dependence of the branching ratio on the values 
of tan (3 in Fig. [| The braching ratios are normalized to the corresponding 
SM values. In the whole parameter range, the branching ratio is greater 
than one, which means that large enhancements can be obtained in the LR- 
SUSY model with respect the SM. This feature is similar to the one in the 
MSSM and is due to the 1/ cos/5 enhancement in the chargino interaction 
vertices. For example, when tan/5 is around 30, an enhancement of one 
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Figure 4: Energy asymmetry A(B — > X s l + l~) normalized to the corresponding SM 
values as a function of tan j3, obtained when m~ g = 300 GeV, fi = 100 GeV, A=50 GeV, 
M L = M R = 500 GeV, m 0q = 300 GeV and m ol = 100 GeV. 

order of magnitude can be obtained. This would make the rare semilep- 
tonic decay more easily to be observed in future experiments. Generally 
the branching ratio increases with tan (3, and for larger values of tan f3 the 
branching ratio will exceed the acceptable range easily. The choice of pa- 
rameters puts stringent restrictions on the allowed values for tan/3. For the 
electron, tan/? should be less than 11 if the branching ratio is to be below 
the experimental bounds, while for the muon tan/? should be less than 26, 
mainly due to the experimental bound of the muon is larger than that of 
the electron. The asymmetries corresponding to different values of tan /3 
are shown in Fig. [|. A clear deviation from the SM model is also obtained. 
Note that in the SM the asymmetry is normalized to 1. The asymmetries 
tend to be large and negative with increasing tan/3. 

We investigate next the dependence of the branching ratio and asym- 
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Figure 5: BK(B — > X s l + l~) normalized to the corresponding SM values as a function 
of m § , obtained when tan/3 = 5, fi = 100 GeV, A=50 GeV, M L = M R = 500 GeV, 
m 0q = 300 GeV and m i = 100 GeV. The experimental constraints are also shown. 
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Figure 6: Energy asymmetry A(B — * X s l + l~) normalized to the corresponding SM values 
as a function of m,g, obtained when tan f3 — 5, n — 100 GeV, A=50 GeV, Ml = Mr = 500 
GeV, m 0q = 300 GeV and m 0/ = 100 GeV. 
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metry on the gluino mass, for a light squark scenario. The chargino and 
neutralino masses are light and fi > 0, a scenario favored by recent anal- 
yses of the anomalous magnetic moment of the muon [|TI| and consistent 
with b — > sj [|l[]. We present the results in Fig. |5]. As the mass of gluino 
increases, the branching ratio will decrease, as the gluino is exchanged as 
a virtual particle in the process. From the branching ratio into electrons, 
the gluino mass is constrained to be heavier than 200 GeV, which is weaker 
than other constraints, for example, from b — > 57; while from the muon 
case there is no constraint on nig. Therefore in the LRSUSY model the 
contributions from gluino-exchanged graphs are not dominant, while this 
is generally so in b — > S7. The corresponding asymmetries are shown in 
Fig. p. It is found that the asymmetries for both the electron and muon are 
less than the corresponding SM value. Although asymmetries do not help 
the experimentalists to observe the decay, they might, if observed, serve to 
distinguish the LRSUSY model from the SM. 

The branching ratio of B — * X s l + l~ is sensitive to the universal scalar 
mass mQ q in the region of small masses only, a feature shared with b — > S7 
[|TJ . This dependence is shown in Fig. [7|. For the electron case, m 0q is found 
to be greater than 200 GeV, where the corresponding scalar quark masses 
are slightly above the current experimental bounds, while for muon case 
there is no constraint on rriQ q . The lepton asymmetries as a function of the 
universal scalar mass, are shown in Fig. ||, and there a small enhancement 
can be found when m 0q is less than 300 GeV. 

In all the previous figures we set the left- and right-handed gaugino 
masses to the same value. In Fig. || and Fig. |10] we investigate the de- 



pendence of the branching ratios and asymmetries on the gaugino mass. 
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Figure 7: BK(B — > X s l + l~) normalized to the corresponding SM values as a function of 
mo, obtained when m § = 300 GeV, tan/3 = 5, fi = 100 GeV, A=50 GeV, M L = M R = 500 
GeV and m i = 100 GeV. The experimental constraints are also shown. 
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Figure 8: Energy asymmetry A(B — > X s l + l~) normalized to the corresponding SM values 
as a function of m , tan (3 — 5, obtained when rrig = 300 GeV, \i = 100 GeV, A=50 GeV, 
M L = M R = 500 GeV and m ol = 100 GeV. 
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Figure 9: BK(B — > X s l + l ) normalized to the corresponding SM values as a function of 
M fl , obtained when m § = 400 GeV, tan /3 = 5, fi = 100 GeV, A=50 GeV, m 0g = 300 
GeV and m / = 100 GeV. The experimental constraints are also shown. 
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Figure 10: Energy asymmetry A(B — > X s l + l ) normalized to the corresponding SM 
values as a function of Mr, obtained when rrig = 400 GeV, tan/5 = 5, fi = 100 GeV, 
A=50 GeV, m 0q = 300 GeV and m ol = 100 GeV. 
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Figure 11: BR(5 — > X s e + e~) normalized to the corresponding SM values versus 
lightest stop mass in the LRSUSY model. 

When Mr = m Qq , the results are not reliable, due to poles in the loop func- 
tions. As the branching ratio of B — > X s l + l~ is dominated by the chargino 
contribution for a large region of the parameter space, one expects a re- 
striction on the left- and right-handed gaugino mass parameters. There 
exist scenarios in which the right-handed symmetry is broken at the same 
scale as supersymmetry, so we expect in those cases to have approximately 



M L = M R Ul9| . With the assumption M L = M R in the gaugino sector, the 
restriction on the right-handed gaugino scale is found to be M R > 400 — 500 
GeV, for low to intermediate squark masses. 



In Fig. 1 I] and Fig. [12] we present the scatter plots of the branching 
ratio and the asymmetry for the decay B — > X s e + e~ as a function of the 
lightest stop mass. We have chosen randomly relevant parameters: tan /3 
changes from 2 to 30, m 0g takes values from 100 to 1000 GeV, \x also varies 
from 100 to 1000 GeV and A = m 0q , while m~ g = 500 GeV, M L = M R = 500 
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Figure 12: Energy asymmetry A(B — ■> X s e + e~) normalized to the corresponding SM 
values versus the lightest stop mass in the LRSUSY model. 

GeV and mo/ = 100 GeV. In addition to the current experimental bounds 
on the SUSY spectra, we also impose the constraints from the rare decay 
b — > 57. It is found that the lightest stop masses lower than 300 GeV are 
excluded. The branching ratio fits easily within the experimental bound for 
large mj, which explains the large number of plot points in that region. An 
enhancement for the branching ratio of one order of magnitude is possible 
while the asymmetry could be 50% larger than the SM value. 



In Fig. [13] we show the correlation between the branching ratio of the 
b — > 57 and B — > X s e + e~ in the above specified parameter ranges. Both 
branching ratios are normalized to the corresponding SM values. Although 
it seems possible that b — > sj is below the SM value, while B —>■ X s e + e~ 
is enhanced with respect to the SM value, there exists a region of the 
parameter space in which both are significantly enhanced with respect to 
their SM values. 
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Figure 13: A correlation between BR(6 — > S7) and BK(B — > X s e + e ) in the LRSUSY 
model. Experimental bounds are also shown. 

4.2 The unconstrained LRSUSY model 

When supersymmetry is softly broken, there is no reason to expect that 
the soft parameters would be flavor blind, or that they would violate flavor 
in the same way as in the SM. Yukawa couplings generally form a matrix 
in the generation space, and the off-diagonal elements will lead naturally 
to flavor changing radiative decays. Neutrino oscillations, in particular, 
indicate strong flavor mixing between the second and third neutrino gener- 
ations, and various analyses have been carried out assuming the same for 
the charged sleptons. In the quark/squark sector, the kaon system strongly 
limits mixings between the first and the second generations; constraints for 
the third generation from b — >• 57 in the LRSUSY model are studied in Ref. 

0- 



We parametrize all the unknown soft breaking parameters, coming mostly 
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from the scalar mass matrices, using the mass insertion approximation (MI) 
f20| . In this framework we choose a basis for fermion and sfermion states in 



which all the couplings of these particles to neutral gauginos are flavor di- 
agonal. Flavor changes in the squark sector arise from the non-diagonality 
of the squark propagators. These off-diagonal squark mass matrix elements 
are assumed to be small and their higher orders can be neglected, and the 
normalized parameters used in the analysis are: 

r ( m d,LL)ij x _ ( m d,RR)ij 
°d,LL,ij — 2 ' °d,RR,ij ~ 2 ' 

m 0q m 0q 



i m d,LR)ij r ( m d,RL)ij , r _, 

°d,LR,ij ~ —2 > O d ,RL,ij ~ —2 , \?() 

" i 0q '"'Oq 



where {fn 2 dAB )ij, A,B = L,R are the off-diagonal elements which mix 
down-squark flavors for both left- and right-handed squarks. We assume 
significant mixings between the second and third generation in the down- 
squarks mass matrix only. We also consider terms with one mass insertion 
only. Although it was shown in the MSSM that double mass insertions 
could possibly enhance the decays of the K meson, their effects on B me- 
son decays are assumed to be negligible. This procedure allows an analysis 
of the graphs contributing to b — ► sl + l~ in terms of a small number of 
parameters. The contribution of each graph with the MI is obtained from 
the constrained case following these simple rules: 

A left gaugino-gaugino vertex has a factor of <^,ll,23 associated with it; a right 
gaugino-gaugino vertex has a factor of Sd^R^ ■ 

A left gaugino Higgsino vertex has a factor of Sd,LR,23 associated with it; a right 
gaugino-Higgsino vertex has a factor of Sd,RL,23- 
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Figure 14: BK(B — > X s e + e~) normalized to the corresponding SM values as a function 
of 5 dtLRi 2 3 , obtained when tan/3 = 5, /i = 100 GeV, A=50 GeV, M L = M R = 500 GeV, 
irtQ q = 500 GeV and moz = 100 GeV. The different lines correspond to different values of 
x = m^/ml q = 0.64(dashed), l(solid), 1.44 (dot — dashed). The experimental constraints 
are also shown. 

• In addition, in the dipole contributions, the term coming from chirality being flipped 
on the fermion leg, proportional to m^/mb or nig/mb, has a factor 8d,LR,23 associ- 
ated with it for the LR contribution, and Sd,RL,23 associated with it for the RL 
contribution. 



With the definition of the mass insertion as in Eq. (p7|), we can inves- 
tigate the effects of intergenerational mixings on the B — > X s l + l~ decays. 
We keep our analysis general, but to show our results, we select only one 
possible source of flavor violation in the squark sector at a time, and assume 



the others vanish. In Fig. [14] we show the dependence of B — > X s e + e~ as 
a function of Sd t iR,23, when this is the only source of flavor violation. The 
horizontal line represents the experimental bound on the branching ratio. 
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Figure 15: BK(B — > X s e + e~) normalized to the corresponding SM values as a function 
of 5 dtLRi 2 3 , obtained when tan/3 = 5, /i = 100 GeV, A=50 GeV, M L = M R = 500 GeV, 
irtQ q = 500 GeV and moz = 100 GeV. The different lines correspond to different values of 
x = m^/ml q = 0.64(dashed), l(solid), 1.44 (dot — dashed). The experimental constraints 
are also shown. 

The branching ratio is plotted as a function of different values for the mass 
ratio x = rn^/m^ q . Fixing rriQ q = 500 GeV, this corresponds to gluino 
masses of 400, 500 and 600 GeV respectively. Constraints on positive and 
negative values of 5d,LR,23 are slightly different, 5±lr,23 is constrained to 
be positive for small mass ratios, and the absolute value of 5d,LR,23 is less 
than 10%. This flavor violating parameter can be strongly constrained 
from b — > 57 because through the Sd,LR,23 term, the helicity flip needed 
for b — *■ sj can be realized in the exchange particle loop. The constraint 
obtained here is complimentary to that from b — *■ sj [0], [21], [22| . The results 
for B — > X s fi + fi~ are not restrictive, so we will not show them here. 

The situation is different when the only source of flavor violation is 



&d,RL,23, as shown in Fig. [15] Again, the most restrictive case is for B 
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X s e + e~ and for the same values of squark and gluino masses as before, only 
positive values of 5d,RL,23 in a small interval close to zero would satisfy the 
experimental bounds, namely 8d,RL,23 < 4%. 

The restrictions on the branching ratio of B — > X s l + l~ from the chirality 
conserving mixings 8d,LL,23 and 8d,RR,2z respectively, with the proviso that 
these are the only off-diagonal matrix elements in the squark mass matrix 
squared, are not as pronounced as the ones for chirality flipping parameters. 
8d,LL,23 and 5^^,23 can almost take all the values in the range (—1.0, 1.0). 
We don't show the results here. 

5 Conclusions 

We analyze the FCNC semileptonic decay B —>■ X s l + l~ in a fully left-right 
super symmetric model. Explicit expressions for all the amplitudes involved 
in the process are given. Constraints on the parameter space of the model 
are obtained in both the constrained case (where the only flavor violation 
comes from the CKM matrix) and the unconstrained case (including soft 
supersymmetry breaking terms). We also include and compare with con- 
straints from b — > S7. 

As a general feature, both b — > 57 and B — > X s l + l~ exhibit similar de- 
pendences on squark and gluino masses. From restrictions on both decays, 
we expect rrig > 250 — 300 GeV and rriq > 200 GeV. A more careful analysis 
of the branching ratio of B — > X s l + l~ reveals that, varying all other pa- 
rameters in the model, the mass of the lightest scalar top should be > 300 



GeV, which is much more restrictive than the experimental bound |23|. The 



parameter that most sensitively affects the branching ratio of B — > X s l + l 
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is tan (3. The constraint from B — > X s / + /~ is slightly more restrictive than 
for b — > S7, and the semileptonic branching ratio can exceed the experi- 
mental bound for t an (3 > 10 — 11, for low squark and gluino masses. In all 
our analysis, we keep \l > 0, and in regions allowed by both (g — 2) M and 
b — ► 57. 

An analysis of the correlation between the branching ratio of B — > X s 7 
and B — > X s l + l~~ reveals that there is a larger region of parameter space in 
which B — > X s / + / _ is enhanced with respect to the SM value by factors of 
almost 10, while B — ► X s 7 is at most 20% larger than the SM value. We 
expect the enhancements to come mostly from regions of intermediate or 
large tan (3. The LRSUSY model shares the strong tan (3 dependence with 
the MSSM, except that here the enhancement is even more pronounced. 
The asymmetry also shows different features from the MSSM ||24j| : it does 



not peak when the Higgs and gauge induced flavor violation are of the 
same size (tan/3 = 35), since the contributions from the gaugino sector are 
different. Also, for low and intermediate values of supersymmetric masses 
(m ,rrig, M L , M R = 200 — 500 GeV), this value of tan/3 is ruled out by 
constraints from B — > X s e + e~ . 

In the unconstrained model, allowing for flavor-dependent soft mixing 
between the second and third generation of squarks (both chirality con- 
serving and chiralty violating), no reliable limits are set on either the LL or 
the RR mixing. However, the chirality violating soft mixing parameters are 
strongly constrained. In particular, the RL mixing, S^rl,23 is constrained 
within four percentage to be close to zero from the bounds on B — > X s e + e~ , 
for a variety of squark and gluino masses; while the constraints on 5d t LR,23 
favor positive values up to 10%. It could be difficult to compare these val- 
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ues with the unconstrained MSSM [p5[| . The bound on Sd,LR,23 appears to 
be much stronger in the unconstrained LRSUSY model than in the MSSM, 
where maximum enhancements are obtained at values of the left-right split- 
ting ruled out in the LRSUSY model. And certainly the restriction on the 
8d,RL,23 coming from B — > X s e + e~ is, to our knowledge, new. 

It appears likely that the most distinguishing factor of the LRSUSY 
model from the SM would be the forward-backward lepton asymmetry. 
This asymmetry, like the branching ratio, is most sensitive to variations 
in tan (3 and could be spectacular even in regions of tan f3 allowed by con- 
straints on the branching ratio. The asymmetry tends to be large and 
negative with increasing tan (3, whereas it is small and positive when vary- 
ing other parameters. As always, the regions of interest are regions of small 
to intermediate values (allowed for branching ratios) for gluino, chargino 
and squark masses. These enhancements are much more pronounced in the 
LRSUSY model than in the MSSM, and increases in the asymmetry by a 
factor of 2 with respect to the SM value are allowed, for a large region of 
the parameter space. 

In conclusion, the decay B — > X s l + l~ would provide an interesting, and 
complementary to h — > sj, test of the LRSUSY model. 
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Appendix 

The relevant Feynman rules and loop functions used in the calculation 
are listed in this Appendix. For further details, we refer to [p]]. The terms 
relevant to the masses of charginos in the Lagrangian are 

f )(:-)+-■• 

where ip + = (-iA£, -iX%, fc, ft, A+) T and V" = (— iA£, -iA£, 0~ 2 , ft, ^) 3 
and 
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where we have taken, for simplification, /Zjj = fi. The chargino mass eigen- 
states Xi are obtained by 

xt = W/, XT = Uij^J, hj = 1, • • • 5, 
with and U unitary matrices satisfying 

WXV' 1 = M D , 

The diagonalizing matrices U* and V are obtained by computing the eigen- 
vectors corresponding to the eigenvalues of X*X and XX*, respectively. 
The terms relevant to the masses of neutralinos in the Lagrangian are 
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where ^0 = ( _ iA 3 ^ _ tX 3 Ri _ tXy ^ ^ ^ ^ ^ ^0 > £b )T and 
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The mass eigenstates are defined by 

X-=N^ (*,j = l,2,...9), 

where iV is a unitary matrix chosen such that 

N*YN~ 1 = N D , 

and N D is a diagonal matrix with non-negative entries. 

In the interaction basis, (q L ,q R ), the squared-mass matrix for a squark 
of flavor / has the following forms. For U-type squarks 

M 2 = ( m 2 + M z (T^-Q u sm 2 6 w )cos2(3 m Uk (A — fi cot (3) \ 



Uk \ m Uk (A- 11 cot 0) 

and for D-type squarks 



rriQ + M Z Q U sin' 1 9 W cos 2f3 I ' 



ml + Mf (Tj - Q d sin 2 6 W ) cos 2(3 m dk (A-p tan f3) 

m dk (A — //tan/3) m 2 , + MfQ^sin 2 W cos 2/? 



The corresponding mass eigenstates are defined as 

_ -pt 

<lL,R — 1 Q L,rQi 



where Tq lr are 6x 3 mixing matrices. The same expressions, with the 
switches Q — * L, U — > iV and D ^ E exist for the sleptons and sneutrinos. 
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The chargino-quark-squark mixing martices G and H are defined as 



GlJL = Vjl{K C KM)il(J'uL)kh 
Gur = Uj 2 (K C KM)il(^UR)kl, 

H 3 UL = 7? H < ^Uj i )(KcKM)ili^UL)kU 

\/2m w sm/3 cos/5 
and the gaugino-sneutrino-lepton G^l,r are defined as 

The neutralino-quark-squark mixing matrices Gq and Hq are defined as 

Gqdl = [sine w Q d N' jl + ^—(T 3 d -Q d sm 2 e w )N' j2 
Vcos 29 w Q u + Q d 

~ Z^Ta O ^j3\y K CKM)il{i- DL)kh 

COS t>w z 

n jki r • n r\ AT' Q d Sin2 ® W AT' 

g odr = -[sm OwQaNj! N j2 

J COS Vw 

+ VC0S *° W (Tj - Q d sm 2 9 w )N' j3 }(K CKM UT DR ) kl , 
cos 9 w 

= 75^ ( S JV - + ^^ )(w) - i(rDl) "' 

and the gaugino-slepton-lepton mixing matrices G el, Gqer are defined 
Gqel = [sm9 w Q e N' jl + ^— {T!-Q e sm 2 e w )N' j2 

J COS (V J 
VCOS 29 W Qe AT , n \ 

Goer = -[sin 9 w Q e N' n - Q ^ ° W N> 2 

J cos 6V J 

+ Vc0S f W (T e 3 - Q e sin 2 9 W )N> 3 ] (T ER ) ki . 
cos 6*^ J 
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The one, two and three variable functions appearing in the decay b — > sl + V 



arc 
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[x 3 + 10a; 2 - 29a; + 18 - (8a; 2 - 6x - 8) log a;], 
[-7x 3 + 8x 2 + llx - 12 - (2a; 3 - 20x 2 + 24x) log a;] 



2(x- 1) 

2 z 

3 x 



- [x 2 - 7x + 6 + (3x + 2) log a;] , 

2 (i + t) V^^lv^}" 1 ' ifx<2 



1 + — 



1 - i In 



i- 



111 



if a; > 2 



a; 



x — 1 



x 



(a; - 1) ; 



x-y 
1 

x-y 
1 

x-y 
1 

x-y 

— I— 

x — y [ x — z 
1 f 1 



log x — (x — > y) 



logx - ^x - (x -»• y) 
x — 1 2 



logx — 



1 



X 



(x-iy 



log x — 



x — 1 

1 



x — 1 



- (x -»• y) 
(x -> y) 



logx x — (x — > z) 

x — 1 2 



x — y l x — z 



X 



x — 1 



logx — (x — ► z) 



- (x ^ y) 
(x -> y)J . 
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